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Abstract - The unbiased thermal diffusion of an overdamped Brownian particle in a square lattice 
potential is considered in the presence of an externally applied ac driving. The resulting diffusion 
matrix exhibits two orthogonal eigenvectors with eigenvalues D\ > D2 > 0, indicating anisotropic 
diffusion along a "fast" and a "slow principal axis". For sufficiently small temperatures, D\ may 
become arbitrarily large and at the same time D2 arbitrarily small. The principal diffusion axis 
can be made to point into (almost) any direction by varying either the driving amplitude or the 
coupling of the particle to the potential, without changing any other property of the system or 
the driving. 



Introduction. — Thermal diffusion plays a key role 
in many physical, chemical, and biological processes. At 
thermal equilibrium, the free diffusion as considered, e.g., 
by Einstein [1] , is generically reduced by an additional pe- 
riodic potential [5] , and - at least for the most common pe- 
riodic lattice structures - remains spatially isotropic (see 
below). Out of equilibrium, the main focus has so far 
been on periodic potentials in one-dimension, perturbed 
by either a static "bias force" (tilted washboard poten- 
tials) [3]-15] or by an unbiased, time-dependent driving [Hj . 
The most remarkable finding in both cases is that the 
diffusion coefficient D may exhibit a giant enhancement 
over the free (Einsteinian) diffusion coefficient Dq. More 
precisely, for asymptotically small Dq, the ratio D/Dq di- 
verges, while D itself still tends to zero, apart from certain 
fine-tuned (non-analytical) potential shapes and system 
parameters, for which D may remain finite. 

Another topic of considerable recent interest concerns 
the effect of disorder (random deviations from a strictly 
periodic potential) , which may give rise to a further diffu- 
sion enhancement or even to anomalous diffusion ^Tj. 

Works on genuine new phenomena in higher dimensions 
are scarce: Guantes and Miret-Artes [8] studied biased dif- 
fusion induced by an asymmetric external driving. Sancho 
et al. and Lacasta et al. [5] put their main focus on (pos- 
sibly transient) anomalous diffusion effects. Experimen- 
tally, Tierno et al. observed giant transversal diffusion of 



paramagnetic particles on an uniaxial garnet film in the 
presence of an oscillating magnetic field [10] . 

Here, we reconsider the prototypical situation of Brow- 
nian motion in a more than one-dimensional periodic po- 
tential and we demonstrate that a simple ac driving gener- 
ically results in anisotropic diffusion, whose direction and 
magnitude depend in a quite intriguing way on the poten- 
tial, the driving, and the particle properties. 

Model. — To keep things as simple as possible, we 
specifically focus on a square lattice in two-dimensions and 
on a square-wave (rectangular) driving. Yet, our main 
results readily carry over to various other drivings and 
more general lattice potentials in two as well as in three 
dimensions. 

Our starting point is the the Langevin dynamics 



AIf{t) + r]f{t) = -VC/(f(t)) + F{t) + ^2rikT£,{t) , (1) 

where M is the mass of the Brownian particle, ry its fric- 
tion coefficient, and r = xSx -\- ycy its position in the x- 
j/-plane, with Cx and &y being the unit vectors in x- and 
j/-dircction. Further, U{r) stands for the spatially periodic 
potential, VU{r) its gradient, and F{t) the temporally pe- 
riodic force. The last term in ([l} accounts for thermal 
noise as usual [2]-[4l|6l[8j|9] : T is the temperature, k Boltz- 



mann's constant, and ^(i) = £,x{t)ex -\- ^yit)ey consists of 
two independent, delta-correlated, Gaussian noises S,xit), 
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Fig. 1: Scliematic illustration of the periodic potential from 
((2]), ([3]) (black contour lines) and of the square- wave driving 
from Q (red double-arrow) which act on the Brownian particle 
according to ([T]) (black dot). The depicted choice of parameters 
L = 1, Q — 1, and a = 35° represents our "standard example" 
in eqs. (O, ^ and in figs.[2]-[Zl One of the periodic potential 
minima is located at the center of the plot. 



As announced, we specifically address the case that the 
periodic potential C/(r) is composed of rotation-symmetric 
on-site potentials t/i(r) on a square lattice with period L, 

U{r)^ Y, Ui{r-[me^ + ney]L) . (2) 

m,n— — oo 

As illustrated by fig. [1] our standard example will be the 
Yukawa potential 



UAf) = Q exp{-|r=l/4L}/|r1 



(3) 



where Q > quantifies the coupling of the particle to the 
potential (e.g. the particle charge in the case of a screened 
electrostatic potential). 

Likewise, the announced square- wave driving is formally 
given by 

F{t) = CaA sign[sin(wt)] (4) 

with amplitude A, frequency w, and direction Ca '■— 
6*2: cos a -f eysina. In other words, it switches between 
A and —A every half-period tt/uj, and a represents the 
angle between the driving force and the x-axis, see fig. [1] 

Yet, as already said, our main qualitative findings turn 
out to be largely independent of the above particular 
choice of on-site potential, lattice structure, and driving. 

To further reduce the number of model parameters, we 
can and will choose the units of length, time, energy, and 
temperature so that 



L = l, uj = \ , r] = 1 



1 



(5) 



and we focus on the simplest and most important case [2]- 
miH] that inertia effects are negligible (ovcrdamped limit). 



r{t) ^ -VU{f{t)) + F{t) + V2T^{t) 



(6) 



We remark that in the presence of an additional, ex- 
ternally applied static bias force, a quite interesting re- 
sponse behavior in the form of a directed net particle mo- 
tion arises [TTKT^. In the following, our main focus will 
be on a quite different issue, namely the unbiased diffu- 
sion of a Brownian particle in the absence of any further 
perturbation. 

DiflFusion matrix. — For symmetry reasons, the 
above specified dynamics rules out any systematic particle 
transport, i.e we are dealing with a purely diffusive, unbi- 
ased Brownian motion, characterized by a 2 x 2 diffusion 
matrix D with matrix elements 



D, 



lim 



{r^{t)rJ{t)) 
2t 



(7) 



where r^, rj G {x, y} (see below ([T])) and where (•) indi- 
cates the ensemble average over many realizations of the 
stochastic dynamics. For ergodicity reasons and due to the 
long-time limit in ([T]), the initial conditions are irrelevant 
and we can without loss of generality focus on 



f(0) = 



(8) 



Since D is symmetric and positive, there are two or- 
thogonal eigenvectors with eigenvalues 



£>! > Da > . 



(9) 



Essentially, a statistical ensemble of (non-interacting) 
Brownian particles, all starting from the origin at time 
t = 0, thus evolves in the course of time into a bigger 
and bigger "particle cloud" of ellipsoidal shape, see fig. 
[2] For large times t, and neglecting local details on the 
scale of the spatial period L, this cloud is quantified by a 
Gaussian probability density [MIS] of variance 2Dit along 
some "fast principal axis" cq := e^ cosQ + Cy sin0, and of 
variance 2D2t along the orthogonal "slow principal axis" . 
The angle Q between "fast direction" and x-axis is thus 
180°-periodic, i.e. 



Q and Q + 180° are equivalent. 



(10) 



Free diffusion (vanishing U{f) and F{t)) according to 
Einstein [1] amounts to Di = D2 = Dq with Dq ~ kT/rj 
for the original dynamics ([1]) and Dq = T for the rescaled 
dynamics ([6|). 

Including the periodic potential U{r), but still without 
oscillating force F{t), the system ([2|)-([6]) is still invariant 
under rotations by 90°, implying isotropic diffusion (Z?i = 
D2). 

Including also the driving F{t), this symmetry is broken, 
implying anisotropic diffusion {Di > D2) in the generic 
case. The remaining two symmetry arguments are: (i) It 
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Fig. 2: Typical examples of diffusive "particle clouds". The 
corresponding (unbiased) Gaussian probability densities are 
completely characterized by the variance 2Dit along the long 
principal axis and by the variance 21)2 1 along the short princi- 
pal axis of the ellipsoidal clouds. In other words, the plot must 
simply be rescaled by \/f in the course of time. Quantitatively, 
every single depicted point represents a numerical solution at 
time t = 10^ of the stochastic dynamics ©-((Gjl with initial 
condition (|8| and parameters Q = 1, q = 35°, T = h ■ 10""", 
and A = 6 (blue), A = 7.84 (green). 



is sufficient to consider a S [0°, 45°]. (ii) If the driving 
direction e^ coincides with a symmetry axis of the square 
lattice (i.e. a = 0° or a = 45°) then ee will be parallel or 
orthogonal to Bq. In any other case (0° < a < 45°) it is 
impossible to predict ee a priori. 

After a couple of unsuccessful attempts one furthermore 
realizes that analytical progress is quite hopeless. There- 
fore we now turn to the presentation of our numerical re- 
sults (next Section), followed by a brief account of the un- 
derlying basic physical mechanisms (Section 'Discussion'). 

Results. — We first exemplify the most interesting 
regime of the parameters Q, A, a and later provide a 
more systematic exploration: fig. [3] depicts the "fast" and 
"slow" diffusion coefficients D\ and D2 together with their 
ratio D1/D2 and the "fast direction" 0. A pronounced 
anisotropy of the diffusion is indicated by the fact that 
Di >• D2. The corresponding "particle cloud" (see fig. 
[2] and below eq. (|9])) thus actually resembles a very nar- 
row "particle beam" . Its orientation may assume any 
value between 0° and (almost) —180° upon variation of 
the driving amplitude A for the specific choice of the re- 
maining parameters used in fig. [3] For other such choices, 
an even larger range of 0-values can be covered. Since 
and -|- 180° are equivalent (see (|T0| ). we can conclude 
that essentially any possible orientation of the "fast" diffu- 
sion direction can be realized upon variation of A, without 
changing any other property of the system or the driving, 
see fig. [21 
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Fig. 3: Eigenvalues Z?i,2 of the diffusion matrix and angle 
Q between first eigenvector and i-axis (see also H9ll0p ) versus 
driving amplitude A, obtained - similarly as in fig. [2] - from 
800 numerical realizations over 10^ time units of the stochas- 
tic dynamics (H))-® with initial condition ((HJ and parameters 
Q = 1, a = 35°, and T = 5- 10"" (blue), T = 5-10"^ (red). In 
spite of the very large final time t = 10^ , the average particle 
displacement along the "slow principal axis" (and possibly also 
along the "fast" one) was less than one spatial period L — 1 
for some yl- values. The corresponding D2 (and Di) values are 
below the range depicted in (b) and (a). In (c), the corre- 
sponding results are marked by dotting the curves. They do 
not represent the "true" limits in ((Tj, since t = 10^ is still too 
small. On the other hand, sufficiently large f-values are way 
beyond our numerical capabilities. In (d) we did not dot the 
corresponding A regions since we are quite confident that the 
depicted results are already very close to the "true" long time 
limit. 



Decreasing the temperature T yields larger D1/D2 ra- 
tios for practically all amplitudes A in fig. [3] for which the 
long time limit in ([7]) could be numerically reached, i.e. 
the diffusion anisotropy gets more pronounced. Whereas 
D2 always decreases upon lowering the temperature, Di 
increases within certain "^-windows" around A = 7.4, 
A = 7.8, A = 8.2 etc., and decreases otherwise. 

According to fig. |4l analogous findings are recov- 
ered upon variation of the potential strength Q from ([3]) . 
Hence, particles with different "charges" Q can be diffu- 
sively "beamed" into different directions very "fast" and 
with very high precision. In other words, a periodic "sur- 
face potential" can act as a very efficient particle sorting 
device. 

A more detailed low-temperature asymptotics oi Di [i = 
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Fig. 4: Same as fig. [3] but for varying coupfing strength Q 
and fixed driving amplitude A = 7. 
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Fig. 5: Eigenvalues Di,2 of the diffusion matrix (O, pre- 
sented as Arrhenius-plot (logDi versus 1/temperature), from 
1000 numerical realizations over 10^ time units of the stochas- 
tic dynamics ((2|-((6| with initial condition ([8]) and parameters 
Q = 1, Q = 35°, and A — 7.77. Dashed: Asymptotic eq. pif) . 




Fig. 6: A pair of transporting attractors, arising due to spon- 
taneous symmetry breaking. Each of the two black trajectories 
represents a solution of the deterministic dynamics ©-((S]) with 
g = 1, yl = 7, a = 35°, and r = 0. The "star" of the lower 
right trajectory indicates the initial position at i = tt and the 
"dot" the final positions at i = 57r, i.e. after 2 temporal peri- 
ods of the square- wave driving (Q. At the end of the depiced 
evolution, the solution has advanced along the y axis by exacly 
two spatial periods of the periodic potential (indicated by the 
grey dots, see also fig. [T]). In other words, we are dealing with 
a periodic attractor with v = 1 and = 90°. Likewise, the 
upper left trajectory advances between i = and t = in into 
the opposite direction. 



1, 2) is provided by fig. [51 evidencing an Arrhcnius-type 
behavior 

A - exp{-E,/T} (11) 

with certain "effective energies" Ei < < i?2- The same 
behavior is recovered within all of the above mentioned 
"^-windows". For all other ^-values, one finds (as ex- 
pected and hence not shown) the same asymptotics as in 
(|TT|) . but now with < Ei < £2- 

We thus can conclude that the diffusive particle mo- 
tion can become arbitrarily "fast" and anisotropic for suf- 
ficiently small temperatures T. Likewise, in comparison 
to the free diffusion coefficient Dq = T, the ratio Di/Dq 
diverges and D2/DQ approaches zero for asymptotically 
small T. 

Discussion. — For an intuitive understanding of what 
is going on, we first focus on the dynamics (HI)-® in the 
deterministic limit (T = 0). Without driving {A = 0) the 
particle readily settles down into one of the local minima 
of the periodic potential ^, ^, see also fig. [T] Upon 
increasing A, all those point attractors evolve via a com- 
plicated sequence of bifurcations into (possibly coexisting) 
periodic, quasi-periodic, or chaotic attractors, whose quite 
intricate details also depend on the values of Q and a 
[TTlfT^ . As might not have been immediately anticipated, 
but in fact is rather plausible at second glance, some of 
the so emerging deterministic attractors give rise to a per- 
manent directed particle motion, sec fig. [51 For symmetry 
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Fig. 7: Deterministic attractors by solving (dJ-® with T — 
numerically for a large number of randomly sampled initial 
conditions r(0). (a) Varying driving angle a and amplitude 
A at fixed coupling strength Q = 1. (b) Varying a and Q at 
fixed A = 7. (c) Varying A and Q at fixed a = 35°. White: 
All numerically detected attractors are non-transporting (van- 
ishing net velocity v = limt^aa r{t)/t). Colored: There is a 
symmetric pair of attractors exhibiting spontaneous transport 
of opposite velocities (cf. fig. [6]). The angle 4> between these 
velocities and the a;-axis is indicated by the different colors. 
The overwhelming majority of the colored pixels corresponds 
to = 90° (red), but also a few blue and green pixels are still 
discernible. In most cases, the modulus of the velocities \v\ was 
found to be of order one (not shown). Insets: Magnification of 
dashed rectangles. 



reasons, such "transporting attractors" always appear in 
pairs with opposite transport velocities by way of spon- 
taneous symmetry breaking. Besides such a pair, there 
may or may not coexist additional (transporting or non- 
transporting) attractors. The actual occurrence of those 
transporting attractors and the concomitant quantitative 
transport velocities v are only accessible by numerical 
means: fig. [7] depicts such numerical results for three rep- 



resentative cross sections through the full, 3-dimensional 
a-A-Q parameter space. 

Turning to small but finite temperatures T, the above 
deterministic attractors become metastable states since 
the thermal noise now induces rare transitions between 
them. For weak to moderate driving amplitudes A, these 
are mainly "hopping" transitions between neighboring po- 
tential wells in x- and y-dircction (cf. fig. [1]) at certain 
rates j^ and jy, respectively. Since the driving F{t) pre- 
dominantly acts along the x-direction for < a < 45°, it 
seems plausible that ^^ > Jy Hence one can infer that 



D2 



L\/2 . 



(12) 
(13) 



Therefore, we can conclude that 9 w 0°. From the 
usual Boltzmann-Arrhenius temperature dependence of 
the rates 'yx,y one finally recovers (fTTj) with < Ei < E2. 
For the rest, the rates jx,y and hence l?i,2 from p2ll3p 
may still depend on A (and likewise on Q and a) in a very 
complicated way. Altogether, these considerations explain 
the main features of fig. 1 up to about A^ 6. 

Next we focus on the "red stripes" in fig. [71 indicating 
two symmetric attractors with spontaneous transport in 
y-direction at certain velocities ±u. Denoting by 7j, the 
rate of noise induced transitions between them, a similar 
calculation as in Ref. [51[T3] yields for the diffusion coeffi- 
cient in 2/-direction the approximation 



(14) 



D, « vy{2^y) , 



whereas in x-direction one obtains similarly as in (|12p the 
approximation 



D, 



Lhx/2 



(15) 



Hence, one can conclude that 6 w ±90°. This provides an 
explanation of eq. (|lip with Ei < < £'2 , an explanation 
of fig. [3] within the above mentioned "A- windows" , and 
an explanation of fig. [S] Moreover, we recognize that 
these "A- windows" are in fact the descendants of the "red 
stripes" in fig. [T] 

So far, we tacitly neglected the possibility that besides 
the transporting attractors there may still coexist further 
non-transporting attractors. Numerically, one finds that 
this is indeed the case within a subset of the colored re- 
gions in fig. [7] (not shown). In such a case, the noise in- 
duced transitions between the various attractors are gov- 
erned by several different rates, resulting in similar but 
more complicated estimates for D12 than in ((T^ - (|15p . 

Analogous generalizations apply for the small parameter 
regions in fig. [7] which are colored differently from red, 
corresponding to a spontaneous deterministic transport 
into a direction substantially different from (j) = 90°. 

Coming back to the remaining A- values in fig. [3] - 
namely those larger than about A = 6 but not contained 
in one of the A-windows - it turns out that, essentially, 
again thermal hopping prevails, but now predominantly in 
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steps of \/2L along the bisectrix y = x, hence 9 « —135°. 
Occasionally, one also encounters more complicated types 
of motion, resulting in even smaller 0-values. 

As expected, with increasing temperature T, all these 
feature become more and more "washed out" in fig. [31 
Especially, this applies to the neighborhood oi A = 6.5, in 
accordance with the tightly nested red and white stripes 
in the corresponding region of fig. [T] 

Analogous considerations apply to fig. |31 Overall, we 
thus can conclude that for sufficiently low temperatures T, 
giant anisotropic diffusion arises within all colored param- 
eter regions in fig. [71 and that pronounced variations of 
the principal diffusion axes are expected whenever crossing 
a border between colored and white regions. 

The entire pertinent parameter regions in fig. [7] are not 
very large as far as their measure is concerned, but still 
very appreciable for instance from the following viewpoint: 
For (almost) any given amplitude A it is possible to find 
angles a and coupling strengths Q for which all those effect 
can be observed [12]. 

Finally, it also seems worth mentioning that in contrast 
to giant diffusion enhancement in one-dimensional peri- 
odic potentials [3J[5], in our present case (i) not only the 
ratio Di/Do, but also Di itself diverges for asymptoti- 
cally small temperatures, and (ii) the pertinent parameter 
regions are of finite measure. 

Conclusions. — Starting with Einstein's ground 
breaking work on Brownian motion [T], the subject of dif- 
fusion has attracted a lot of attention due to its practical 
relevance in numerous specific systems, but also due to 
its theoretical interest as a fundamental form of trans- 
port per se. Here, we reconsidered diffusion in a partic- 
ularly simple "minimal model": a spatially periodic sys- 
tem, driven out of equilibrium in one of the experimen- 
tally most natural and common ways, namely by an os- 
cillating driving force. Our first main finding is an ex- 
treme anisotropy of the diffusion, whose "speed" in one 
direction may grow exponentially fast as temperature de- 
creases, while the perpendicular "speed" approaches zero 
exponentially fast. Our second main finding is that the 
principal axes of the anisotropic diffusion can be made to 
point into (practically) any direction by solely varying e.g. 
the ac driving amplitude or the coupling of the particle to 
the periodic potential, but without changing any other 
property of the system or the driving. 

The indispensable prerequisites for those diffusion ef- 
fects are very weak and ubiquitous: Brownian motion in 
a spatially periodic lattice potential under the action of 
an oscillating driving force is a paradigmatic scenario in 
a large variety of different contexts, e.g. diffusion on the 
surface [8] but also in the bulk of a crystal under the action 
of electromagnetic perturbations, migration in periodic ar- 
rays of laser traps, or in periodically structured micro- or 
nano-fluidic systems, to name but a few [51 [TII1[T^ . 

One reasons why those very fundamental diffusion phe- 
nomena have not been discovered earlier may be the fact 



that the underlying key mechanisms of spontaneous sym- 
metry breaking has not been naturally associated with 
such systems up to now. A second reason is that a very 
careful numerical exploration is required with a clear idea 
of what one is actually looking for. In turn, once the quan- 
titative numerical predictions for a specific experimental 
system are available, verifying and exploiting the predicted 
effects seems quite feasible and worth while to us. 
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